Interplay between strong correlations and magnetic field in the 
symmetric periodic Anderson model 



Debabrata Parihari and N. S. Vidhyadhiraja 

Theoretical Sciences Unit, 
Jawaharlal Nehru Centre for Advanced Scientific Research, 
Jakkur, Bangalore 560064, Indiau 

David E. Logan 

Oxford University, Chemistry Department, 
Physical and Theoretical Chemistry Laboratory, 
South Parks Road, Oxford 0X1 3QZ, UK. 
(Dated: June 10, 2008) 

Magnetic field effects in Kondo insulators are studied theoretically, using a local moment approach to the 
periodic Anderson model within the framework of dynamical mean-field theory Our main focus is on field- 
induced changes in single-particle dynamics and the associated hybridization gap in the density of states. 
Particular emphasis is given to the strongly correlated regime, where dynamics are found to exhibit universal 
scaling in terms of a field-dependent low energy coherence scale. Although the bare applied field is globally 
uniform, the effective fields experienced by the conduction electrons and the /-electrons differ because of 
correlation effects. A continuous insulator-metal transition is found to occur on increasing the applied field, 
closure of the hybridization gap reflecting competition between Zeeman splitting and screening of the f- 
electron local moments. For intermediate interaction strengths the hybridization gap depends non-linearly 
on the applied field, while in strong coupling its field dependence is found to be linear. For the classic Kondo 
insulator YbBi2, good agreement is found upon direct comparison of the field evolution of the experimental 
transport gap with the theoretical hybridization gap in the density of states. 
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INTRODUCTION 



Kondo insulator materials such as SmB6, YbBi2 and 
Ce3Bi4Pt3 have been of sustained interest to experimen- 
talists and theorists for several decadesi^^^i Inter- 
est in these system stems from their unusual properties, 
such as the small hybridization gap (of a few meV) and 
mixed valency, as well as the transition to metallic behav- 
ior with doping 8 , application of pressure 9 and magnetic 
field 10 . The underlying qualitative origin of such rich 
behavior is the strong electronic correlation arising due 
to the localized and narrow /-orbitals of the rare earth 
atoms, which hybridize weakly with broad, essentially 
non-interacting conduction bands. 

A quantitative understanding of the dynamics and 
transport properties of these materials has not however 
been easy to achieve. A principal stumbling block in this 
regard has been the theoretical treatment of localized 
and itinerant fermionic degrees of freedom on a compa- 
rable footing. Much progress has been made in recent 
years with the advent of dynamical mean field theory 
(DMFT) n i 12 i 13 i 14 ; within which generic lattice-fermion 
models such as the Hubbard or the periodic Anderson 
model have found approximate solutions, and quanti- 
tative agreement with experiments has also been ob- 
tainedii^. Within DMFTi2*i2*12*ii, a lattice fermion 
model is mapped onto an effective single-site correlated 
impurity which hybridizes with a self-consistent conduc- 
tion electron bath. Thus, various techniques such as the 
numerical renormalization group, exact diagonalization, 



diagrammatic perturbation theory based approaches, 
quantum Monte Carlo and many others that have, in the 
past, been developed to handle the many-body single- 
impurity problem, have now been adapted and modi- 
fied for use within the DMFT framework 1 UU&li. One 
such recent technique is the local moment approach 
(LMA)^^^^^^^^, which has been shown to be 
powerful not only in the context of single-impurity sys- 
tem a 16 i 17 i 18 i 19 , but also for lattice-based heavy fermion 



systems 



15,20,21,22 



when used in conjunction with DMFT. 



In this paper, we employ the LMA to understand the 
interplay of electronic correlations and an external mag- 
netic field in Kondo insulator materials. 

The generic model used to study Kondo insulator ma- 
terials is the periodic Anderson model (PAM), which 
consists in physical terms of a correlated /-level in each 
unit cell hybridizing locally with a non-interacting con- 
duction band. Magnetic field effects in these systems 
have been studied theoretically through the inclusion of 
a Zeeman term in the PA M 23 i 24 . The observed insulator- 
metal transition^ has also been reproduced qualitatively 
in theoretical calculations^. However, a detailed un- 
derstanding of the changes in single-particle dynamics 
and the associated hybridization gap has on the whole 
been lacking, and a quantitative description of the ex- 
perimentally observed field-induced behavior has not 
been achieved. We seek to bridge these gaps in this 
paper by studying the PAM with a Zeeman term us- 
ing LMA+DMFT, and with particular emphasis on the 
strongly correlated (or strong coupling) regime. Our pri- 
mary focus is on the field-induced changes in the single- 



2 



particle dynamics and the associated hybridization gap 
in the density of states. 

The outline of the paper is as follows: We begin in the 
next section with a brief description of the model (PAM), 
the DMFT framework and the LMA technique for the 
PAM in the presence of a magnetic field. In section 3 
we present our theoretical results and their analysis. We 
also make comparison between theory and experiments 
on the classic Kondo insulator material YbB 12. Brief 
conclusions are given in section 4. 



II. MODEL AND FORMALISM 

The Hamiltonian for the PAM is given in standard 
notation by 



+V J2(flc ia + h.c) + e c £ c\ a c ia (1) 



where the first term describes the kinetic energy of the 
noninteracting conduction (c) band due to nearest neigh- 
bour hopping tij . The second term refers to the /-levels 
with site energies €f and on-site repulsion J7, while the 
third term describes the c/ f hybridization via the lo- 
cal matrix element V. The final term represents the c- 
electron orbital energy. Within DMFT— which 
is exact in the limit of infinite dimensions, the hopping 
term is scaled as oc t*j\fZ~ c , with coordination num- 
ber Z c — » oo.In this paper we consider mainly the hy- 
per cubic lattice, for which the non-interacting density of 
states is a Gaussian 11 (po(e) = exp(— e 2 /t 2 )/ y/nt*). We 
also consider specifically the (particle-hole) symmetric 
PAM, which is the traditional limit employed to study 
Kondo insulators^^^^^^^L 2 ^. For the sym- 
metric PAM the conduction band is located symmetri- 
cally about the Fermi level (i.e. e c = 0), while e/ = 
— U/2. This corresponds to half- filling of the / and c- 

levels, i.e. n f = Y Mia ha) = 1 and n c = E<r( c LO = 
1 for all U. In this case the system is an insulator for all 
interaction strengths U (in the absence of a magnetic 
field), with a corresponding gap in the single-particle 
spectra. The symmetric PAM in the absence of a mag- 
netic field has been studied quite comprehensively within 
the DMFT framework, see e.g.U^^ 2 ^ 2 ^. 

Within DMFT, the lattice fermion model maps onto an 
effective, correlated single-site impurity hybridizing self- 
consistently with a conduction electron bath» i 12 i 13 i 14 . 
The self energy is thus spatially local, i.e. momen- 
tum independent. Thus, the problem is simplified to a 
great extent. Nevertheless, the problem remains non- 
trivial because the impurity model is as yet unsolved 
for an arbitrary hybridization. As mentioned in the 
introduction, the local moment approach^ has been 
successful in describing the single impurity Anderson 
mode l 16 i 17 i 18 i 19 as well as in understanding the PAM 



within DMFT— »2Sl»2i*22. Here we extend the approach to 
encompass finite magnetic fields in the symmetric PAM, 
enabling study of magnetic field effects in Kondo insula- 
tors. 

The local moment approach begins with the symme- 
try broken mean field state (unrestricted Hartree Fock, 
UHF). Dynamical self energy effects are then built in 
through the inclusion of transverse spin fluctuations for 
each of the symmetry broken states. The most important 
idea underlying the LMA at T = is that of symmetry 
restoratio n 16 ! 17 ! 22 : self-consistent restoration of the bro- 
ken symmetry inherent at pure mean field level, arising 
in physical terms via dynamical tunneling between the 
locally degenerate mean- field states, and in consequence 
ensuring correct recovery of the local Fermi liquid be- 
havior that reflects adiabatic continuity in U to the non- 
interacting limit. The reader is referred to our earlier 
paper s 15 ! 22 for full details of the formalism and imple- 
mentation of the zero-field LMA for the PAM. 

In the presence of a global magnetic field, the degen- 
eracy between the mean field solutions 16 (labeled as A 
and B corresponding to and - found at the mean 
field level is lifted, and only one solution as determined by 
sgn(h) remains^. Here, we consider explicitly ft > for 
which only the 'A' solution survives. For ft ^ the bare 
electronic energy levels, e a for a = c and /, are of course 
split via the Zeeman effect. The modified energy levels 
are given by e aa = e a — ah a with h a = ^g a fiBH (and 
a = ± for j / [ spins), where \±b is the Bohr magneton, g a 
the Lande g- factor and B the magnetic field. Although 
gf 7^ g c in general, for simplicity we set gf = g c = g i.e. 
we consider the application of a uniform magnetic field 
ft = ft c = hf to both c and /-levels. The local (site- 
diagonal) Green functions for the c- and /-levels may be 
expressed using the Feenberg renormalized perturbation 
theory^- 31 as 



G c Aa (u;h) = [cj + + aft 



V 2 



+ ah — YjfAcr(w, ft) 

SA^G^h)}]- 1 (2) 



G f Aa (oj;h) = [w+ + ah-t fAa (Lu;h) 



V 2 



>+ + ah-S Aa [G c Aa (u J] h)] 



-\ (3) 



where uo + = uo + isgn(u)0+ . Here Sao- is the Feenberg 
self energy, a functional solely of G c Acr1 given by 



S A a(u;h) = jAa(u;h) 



with 



7Acr(^5 ft) = + Crft ■ 



H[ 7 Aa(^;ft)] 



V 2 



uj+ + ah — YifA(j(w] ft) 
and the Hilbert transform H[zl defined as 



m = / 



po(e)de 



(4) 



(5) 



(6) 
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FIG. 1: Principal contribution to the LMA {E f(T (uj)}. Wavy 
lines denote U. See text for details. 

such that G c Aa (uo] ft) = H[jAa]- 

Dropping the 'A' subscript for clarity, the spin- 
summed Green functions are denoted by G a (co]h) = 
\ J2 a G%(uj; ft) (a = c, /), with corresponding spectral 
functions D a (uj;h) — — ^sgn(uj)ImG a (uj;h). Within the 
LMA the /-electron self energies Yif a (u; ft) are familiarly 
separated into a static mean-field contribution plus a dy- 
namical part H/ct(cj; ft)A5*ifi*22, 

Z fa (u-,h) = -lu\fi(h)\ + ZML>-,h) (7) 

where \fi(h)\ is the UHF local moment. We approximate 
the dynamical part of the self energy by the usual non- 
perturbative class of diagrams retained in practice by the 
LMA (as shown in figure [1]), which may be expressed 
mathematically at zero temperatur o 15 i 16 i 22 as 

£ /c >; h) = U 2 [ g- a (u - ftjn-^^i; ft). 

(8) 

Here the host/medium Green function Q G (denoted by 
the double-lined propagator in figure [1]) is defined by 

g a (co; ft)" 1 = GUu; ft)" 1 + E /ff (o;; ft) . (9) 

U~ acr (uj; ft) denotes the transverse spin polarization prop- 
agator (shown shaded in figure [T]), which in the random 
phase approximation employed is expressed as II" 1 = 
°n+-/(l - f7°n + -). The bare polarization propagator 
°n +_ (cc;; ft) is expressed in terms of mean-field propaga- 
tor o 15 i 22 {g%(cu;h)} (in which only the static Fock con- 
tribution to the self-energies occurs). 

For ft = the symmetry restoration (SR) condition for 
the symmetric PAM is given b y 15 i 22 S/ a (cc; = 0; ft = 0) = 
(independent of spin a, since by particle- hole symmetry 
t fa (u; ft) = -t f - a (-uj; ft)), i.e. by 

E /t (cj = 0; ft = 0) = U\p,\/2. (10) 

Satisfaction of the SR condition ensures adiabatic con- 
tinuity (in U) to the non-interacting limi t 15 i 22 of the 
hybridization-gap insulator (the system being in that 
sense a generalized Fermi liquid 15,22 ); and as expected 
the system remains insulating, with a gap in the single- 
particle spectra, for all interactions U > 0. 



The practical implementation of the above procedure 
is carried out as follow o 15 i 21 i 22 , in which the numerical 
procedure is simplified by keeping x = \U\fi\ fixed, while 
varying U to satisfy the symmetry restoration condition. 
We begin with an ft = calculation for a given x. (i) 
The mean- field Green functions g^(uo;h) are obtained 
from equations (2-7) by retaining only the static part 
of the self energy, (ii) These Green functions are then 
used to construct the bare polarization bubble °II H . (hi) 
The DMFT iterative procedure starts with a specific dy- 
namical self energy, obtained either as a guess or from 
the previous iteration, which is substituted into equa- 
tions (2-7) to get the c- and /- Green functions, (iv) 
The host /medium Green function Q a (uj;h) is obtained 
through equation (9). (v) The transverse spin polariza- 
tion propagator II -1 is computed for a given U and along 
with Q a (uj;h), substituted in equation (8) to determine 
the ou = self energy for a given U; (vi) The SR con- 
dition (equation (10)) is checked, and if found not to be 
satisfied, U is varied and step (v) is repeated until the 
SR condition is satisfied, (vii) Upon restoration of the 
broken spin symmetry, the corresponding U is obtained 
and equation (8) may then be used to get the full dy- 
namical self energy at all frequencies, (viii) The new self 
energy is fedback in the first step of the DMFT proce- 
dure (step (hi)) and the iterations are continued until full 
self-consistency is achieved. For finite fields the same pro- 
cedure is adopted except that SR is no longer imposed 18 
(the U found from SR at ft = is naturally used for all 
ft>0). 

III. RESULTS AND DISCUSSION 

Before discussing the interplay between interactions 
and magnetic fields in the symmetric PAM, we consider 
briefly two limiting cases, (i) First the non-interacting 
limit, U = 0, with ft / 0; and secondly (ii) the inter- 
acting problem U > in the absence of a field. In the 
former case, the non-interacting c-spectrum for spin-a is 
given by 

d^;h) = P0 ^ + ah--^j . 

The spectral band edges are given by 

V 2 

uo + crh - = ±W (11) 

u + ah y J 

where 2W is the full band- width of the non- inter acting 
density of states po(cu). It is easy to see from this that 
for ft = there is a hybridization gap at the Fermi level 
(denoted by Ao(0)), which decreases with increasing field 
ft and eventually closes at a field value that is half of the 
zero field spectral gap, i.e. 

Ao(ft) = A (0)-2ft. (12) 

Thus in the non-interacting limit, application of a mag- 
netic field leads to an insulator-metal transition simply 
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FIG. 2: LMA conduction electron spectra for h = 0, in strong 
coupling (left panel, U/t* = 7.1) and weaker coupling (right 
panel, U/t* — 1.2), with V 2 — 0.2£*. Insets show the low 
frequency region of the respective main panels. 



because of the rigid crossing of the up- and down-spin 
bands. We add here that the non-interacting Gaussian 
density of states characteristic of the hypercubic lattice 
(as considered explicitly below) is of course unbounded 
and as such does not possess 'hard' band edges. The field- 
induced insulator-metal transition in this case is thus 
strictly a crossover, although in practice the transition 
is 'sharp' as one would expect (see e.g. figure [5]). 

In the second limit, of finite interactions but zero 
field, the ground state remains gapped for all interaction 
strengths 2 ^, the hybridization gap A decreasing contin- 
uously with increasing U from its non-interacting limit 
Ao(0). In the strong coupling, Kondo lattice regime of 
the model, universal scaling occur o 15 i 22 i 26 i 28 i 29 in terms of 
an exponentially small low energy scale ul = ZV 2 /t*(= 
~A);Z being the quasiparticle weight or mass renormal- 
ization factor, given by Z = (1 — dY i f a (u)/du\ u =o)~ 1 . 
The Green functions and their associated spectra depend 
solely on u = u/ul in the universal scaling regim o 15 i 22 . 
Representative results for the zero-field density of states 
are shown in figure [2j where the main panels show the 
conduction band density of states, t*D c (u) (solid lines), 
as a function of frequency, u/t*. The right panels rep- 
resent intermediate coupling (U/t* = 1.2), while the left 
panels are for strong coupling {U/t* = 7.1). The insets 
show a close-up of the low frequency spectra where the 
hybridization gap at the Fermi level (u = 0) is evident 
at both weak and strong coupling. 

Now we consider both interactions and the field. In 
the strong coupling regime, in parallel to the h = 
limit, we expect universality to persist in terms of a 
field-dependent low energy scale, ul = ul(Ji). To de- 
rive explicitly the universal scaling form in the limit of 
low-frequencies (i.e. close to the Fermi level), we perform 
a simple low-frequency 'quasiparticle expansion' of the 
self energy, retaining only its real part (^f a ) to leading 
order in u: i.e. 



EUu-h) = Ef a (0-h) 



1 



Z(h) 



1 u 



(13) 



where Z(h) = (1 - dY.f a (u; h) / du^y 1 is the field- 
dependent quasiparticle weight (independent of a since 



T,^ a (u; h) — —T^_ a [—u\ h) by particle- hole symmetry). 

Substituting equation (13) into equations (0|7j), we find 
that the associated spectral functions are just renormal- 
ized versions of their non-interacting counterparts, being 
given by 



D c a (u;h) ^-4° po 



u + crh e ff 



Dt(u;h) 



u>^o ti 1 



V 2 (<5 + Ghefff 



(14) 



D c M,h) (15) 



where u = u/ul(K), and the low-energy scale ul(K) = 
Z(h)V 2 jt* is thus defined (in direct parallel to the h = 
limit). In obtaining equations (|14|15p we have explicitly 
considered the strong coupling scaling regime, of finite 
u = u/ulQi) and h/uL(h) in the formal limit where the 
low-energy scale ul — > (so that 'bare' factors of u = 
uul or h are thus neglected). h e ff in equations (|14I15|) is 
given by 



eff 



h 



V 2 



(16) 



(being independent of a, by symmetry); or equivalently, 
using the symmetry restoration condition 1^(0; 0) = 0, 

by 



h + a(t%(0;0)-tf a (0;h)) 



V 2 ' 



(17) 



In physical terms, h e ff represents a dimensionless effective 
field, and its primary field-dependence arises from that 
of the interaction self-energy (the 'bare' factor of h in 
equations ([T6]) or (|T7|) can of course be dropped in the 
strict scaling limit, although we retain it for clarity). In 
fact a leading order Taylor expansion of equation ([T6]) 
or (fT7| gives h e ff = c/i, where h = h/uL(0) is the field 
rescaled in terms of the h = low energy scale ul(Ji = 
0) = Z(0)V 2 /t*, and c = Z(0)/Z (~ (9(1)) with Z = [1- 
a(dT l ^ a (0]h)/dh)h=o}~ 1 thus defined. From this simple 
consideration we anticipate that h e ff is just a rescaled 
version of the bare magnetic field itself, and is on the 
order of h = H/ul{^) (as confirmed explicitly below, see 
figure [5]) . 

Equations (jEj) and (fT5|) show that in strong coupling, 
the spectra D c (u;h) and V 2 D^ (u^h) should be univer- 
sal functions of a; = u/uL(h) : for a fixed h e jf. Thus, 
if distinct sets of model parameters in the strong cou- 
pling regime correspond to the same h e ff , the spectra D c 
and V 2 Df should collapse to the same scaling form as a 
function of u, independently of the bare parameters U/t* 
and V/t*. That this is is so is illustrated in figure O 
where the top panel shows the full LMA c-electron spec- 
tra t*D c for the hypercubic lattice, and the bottom panel 
the corresponding /-electron specta (V 2 /t*)D^ . Three 
sets of spectra are shown, with parameters Ut*/V 2 = 35 
and h/t* = 2.4 x 10~ 5 (dashed), Ut*/V 2 = 30 and 
h/t* = 7.4 x 1(T 5 (dotted) and Ut*/V 2 = 25 with 
h/t* = 2.3 x 10 _4 (solid); in each case, h e jf (= 0.27) is 
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FIG. 3: LMA spectra for the c-electrons (top panel, 
t*D c (uj; h)) and /-electrons (bottom panel, V 2 t~ 1 D^ (cj; h)), 
for three parameter sets, each corresponding to a fixed h e ff = 
0.27: UU/V 2 = 35,/i/t* = 2.4 x 10" 5 (dashed), UU/V 2 = 
30, h/U = 7.4 x 10" 5 (dotted) and UU/V 2 = 25,/i/t* = 
2.3 x 10 -4 (solid). The insets show the spectra vs. 'bare' 
frequency, u/t*. Main panel: The same spectra plotted vs. 
Q = u>/u>L(h) collapse to a common scaling form. 



the same. The insets to the figure show that the spectra 
as a function of the 'bare' frequency uj/t* are distinct. 
However when plotted vs. Co (as shown in the main pan- 
els), they are indeed seen to collapse to a single universal 
form. 



The quasiparticle forms in equations (|14|15p embody 
local Fermi liquid behavior and adiabatic continuity to 
the non-interacting limit. They give explicitly the leading 
low-frequency asymptotic behavior of the scaling spec- 
tra that must be satisfied by any 'full' theory. Direct 
comparison between the quasiparticle forms and the full 
LMA scaling spectra is shown in figure[l]for the c-electron 
spectra, and for two values of the effective field h e ff (one 
corresponding to a case where the system remains in- 
sulating, the other for a higher field where the system 
is metallic, as discussed below). It is indeed clear from 
the figure that the LMA correctly recovers the limiting 
quasiparticle form in the vicinity of the Fermi level. In 
physical terms it is also worth noting that the low fre- 
quency quasiparticle spectra are essentially those for the 
non-interacting limit, but with the local fields for /- and 
c-electrons replaced by (V 2 /t*)h e ff and zero respectively. 
Hence, although the bare applied field is globally uni- 
form, the effective local fields experienced by the c and 
/ electrons are different because of correlation effects. 

We turn now to the transition with increasing field 
from an insulating state characterised by a spectral gap 
straddling the Fermi level, to a metal with a finite density 
of states at uo = 0. Equations PHI and [T5l may be used to 
obtain an estimate of the spectral band edges in strong 
coupling, and hence the gap as a function of the field. 




FIG. 4: Full LMA c-electron scaling spectra (solid lines) are 
compared to the limiting quasiparticle form equation ([14)) 
(dashed lines). The left panel shows the scaling spectra for 
h e ff =0.27 (where the system is an insulator) and the right 
panel for h e jf — 0.51 (metallic). 
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FIG. 5: c-band spectrum at the Fermi level D c (0;h) in the 
strong coupling, universal regime as a function of h e ff, showing 
a continuous insulator-metal transition at h e ff, c — 0.36. Top 
inset: shows the linear dependence of the effective field h e jf 
onh = h/oo L (0) (vl(0) = Z(0)V 2 /U). Bottom inset: linear 
increase of the quasiparticle weight Z(h) with increasing h. 



The band edges are given by 
1 



Co + crh e ff 



= ±w/u 



(18) 



with 2W the band-width of the non-interacting spec- 
trum. From this the field-dependent gap in D c (uj;h) or 
Df(u;h) follows as 

i( ^ 2(l _^)»f. (19) 

This in turn implies an insulator to metal transition at a 
critical effective field h e ff iC that is on the order of unity 
(~ t*/W). The main panel in figure [5] shows the vari- 
ation of the full LMA density of states at the Fermi 
level, D c (ou = 0; h) (calculated explicitly for U/t* = 
6.1,F 2 /tJ = 0.2), as a function of h e jf on a log scale. 
Although the calculations are for a hypercubic lattice, 
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with a strictly soft gap in its zero-field spectrum, the 
insulator-metal transition is seen in practice to be sharp; 
occurring at a critical ft e jf, c ~ 0.36 that is indeed on the 
order of unity (we identify the critical field in practice 
from UD c (0;h c ) ~ 10~ 3 ). On further increase of the 
field, D c (0; ft) is seen from the figure to rise continuously, 
towards the high- field value of l/y^rt* which is just the 
non-interacting dos value at the Fermi level. 

The top inset to figure [5] shows the dependence of the 
effective field h e jf (equations (|16|17p ) on the scaled ex- 
ternal field ft = h/cu L (0) (with u L (0) = Z(0)V 2 /U). h eff 
is seen to be linear in ft (which behavior extends over a 
wide ft interval) and, as anticipated above, is of the same 
order as it: h e jf ~ ft/2 as evident from the figure. The 
lower inset to the figure also shows the /^-dependence of 
the quasiparticle weight Z(h). It too is seen to increase 
linearly with field, implying a lowering of effective mass 
with an increase in field; and which behavior is consistent 
with a similar finding for the single impurity Anderson 
models. 

The field-dependence of the full density of states is il- 
lustrated in figure [6l where we plot the (universal, strong 
coupling) conduction band density of states D c (u; ft), as 
a function of the &, for various h e ff. The solid curve 
h e ff = represents the insulating ground state, while the 
dotted curve is for h e ff = 0.38, which is just above the 
insulator-metal transition, so the gap has closed. The 
remaining curves are for h e ff = 1 and h e ff = 5, show- 
ing metallic densities of states characterised by a finite 
spectral density D c (0; h) at the Fermi level. 

0.8 1 i i i ■ i ' i ' i I 
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FIG. 6: Universal c-electron spectra from the LMA are shown 
as a function of Co — oj/ul for various fields: h e ff = (solid), 
heff = 0.38 (dotted), h e ff = 1 (dashed) and h e ff = 5 (dot- 
dashed). The closure of the insulating gap with increasing 
field is evident. 

In the non- interacting limit, [7 = 0, the spectral gap 
closes linearly with the applied field as in equation [T2l 
and the essential mechanism for the insulator-metal tran- 
sition is obvious: Zeeman splitting moves the up- and 
down-spin bands rigidly, resulting in their crossing at a 



critical field, ft c o/Ao(0) = \. This simple picture is nat- 
urally modified in the presence of correlations, U > 0, 
where two essentially competing effects are operative. 
First, the tendency of the system to lower its energy by 
uniform ('ferromagnetic') spin polarization of the c- and 
/-electrons, i.e. the Zeeman effect, which alone operates 
in the non-interacting limit. However for h = in the 
presence of interactions, lattice-coherent Kondo singlet 
formation occurs, driven by local anti ferromagnetic spin 
correlations between the c- and /-electrons. In the pres- 
ence of both interactions and a field, Zeeman splitting 
thus in effect competes with local moment screening; and 
the field-dependence of the spectral gap is not a priori 
obvious. 




FIG. 7: The field-dependent spectral gap scaled by the zero 
field gap, A(ft)/A(0), vs. h/h c , for U/U = 1.2, V 2 /^ = 0.2 
(triangles), U/U = 5.1,V 2 /tl = 0.2 (squares) and U/U = 

6.1, V 2 /tl = 0.2 (circles). (The lines are best fits to the 
points.) Inset: the gap in units of t* vs. the bare field h/t*. 

LMA results for the spectral gap are shown in figure [71 
where the field dependent gap scaled by the zero field 
gap, A(ft)/A(0), is plotted vs. h/h c for various inter- 
action strengths. For intermediate coupling (U/t* = 

1.2, V 2 /t 2 = 0.2 (triangles)), the gap is seen to close 
non-linearly in the field and is best fit by a quadratic 
form. In the strong coupling regime by contrast (squares, 
U/U = 5.1, V 2 /tl = 0.2 and circles, U/U = 6.1, V 2 jt\ = 
0.2), linear behavior A (ft) = A(0)(1 — h/h c ) is obtained, 
similar to the non-interacting limit. In this case how- 
ever, when A(ft) is plotted directly vs. the bare field 
h/t* as shown in the inset of figure [3 the functional 
form obtained is A (ft) = A(0) — ft; showing that the 
field ft c required to close the gap in strong coupling sat- 
isfies ft c /A(0) = 1, i.e. twice that required in the non- 
interacting limit, where ft c o/Ao(0) = |. This result is 
physically natural, in view of the effective competition 
between Zeeman splitting and local moment screening 
discussed above. 

Finally, we would like to make a comparison of our 
theoretical results to experiment. For the classic Kondo 
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insulator YbBi 2 , the field dependence of the transport 
gap has been determined from low-temperature resistiv- 
ity measurements 10 (the leading low-T behavior of the 
resistivity being p(T) oc exp(— A tr /T) with A tr = A tr (h) 
the transport /activation gap). Since the transport gap 
is known theoretically 15 to be proportional to the spec- 
tral gap (A(0) — 2A £r (0)^ 5 -), comparison to experiment 
may be made. In our earlier work^ 5 - where we compared 

I ■ i • i 1 i I 




q| , 1 , j , 1 

150 300 450 

h [kOe] 

FIG. 8: Comparison of the experimental 10 transport gap in 
YbBi2 (opencircles) to the theoretical gap (filled squares), 
obtained as discussed in text. The theory is seen to describe 
well the functional form of the transport gap. 

zero-field transport properties of YbBi2 to theoretical re- 
sults from the LMA, we concluded that YbBi2 belongs to 
the intermediate coupling regime (and as such lies out- 
side the universal scaling regime). This is corroborated 
by the field-dependence of the transport gap, experimen- 
tal results for which 10 are shown as open circles in fig- 
ure [8j The dependence of A tr (h) on the field h is clearly 
non-linear, which behavior we have found above to be 
characteristic of the intermediate coupling regime. To 
make comparison to experiment in this regime, specific 
model parameters must of course be specified, and here 
we choose U/t* = 1.2, V 2 /t 2 = 0.2 (the essential results 
are quite insensitive to these particular values). The filled 
squares in figure [8] show the field dependence of the resul- 
tant theoretical spectral gap, compared directly to exper- 
iment with a simple multiplicative scaling of the x and y 
axes. The functional form of the theoretical gap is seen 



to be almost identical to that found experimentally, thus 
yielding good agreement between theory and experiment. 
Further, since the experimental A tr (0) ~ 60K (figured]) 
then the spectral gap A(0) ~ 2A £r (0) ~ 120K; and for 
the bare parameters considered we find A(0) = 0.026£*. 
This in turn yields the estimate t* ~ OAeV, which is 
physically realistic and compatible with transfer integral 
values found through a band structure calculation 32 . 
IV. CONCLUSION 

The interplay between electronic correlations and an 
externally applied magnetic field in Kondo insulators 
has been considered in this paper. The symmetric pe- 
riodic Anderson model, with a Zeeman term to account 
for the external magnetic field, has been studied within 
the dynamical mean field framework using a local mo- 
ment approach. In the strong coupling Kondo lattice 
regime of the model, the local c- and /-electron spec- 
tral functions are found to exhibit universal scaling, be- 
ing functions solely of uo/uol (with ujL(h) = Z(h)V 2 /t* 
the characteristic low-energy scale) for a given effective 
field h e ff. Although the externally applied field is glob- 
ally uniform, the effective local field experienced by the 
c- and /-electrons differs because of correlation effects. 
The zero-field spectral gap characteristic of Kondo insu- 
lators is found to close continuously, leading to a contin- 
uous insulator-metal transition at a critical applied field 
h c . Field induced closure of the insulating gap is not 
simply a rigid band-crossing affair, but involves competi- 
tion between local moment screening (reflecting correla- 
tion effects) and Zeeman spin-polarization. In the inter- 
mediate coupling regime the gap is found to close non- 
linearly with field, while in the strong coupling regime it 
closes linearly. Comparison of the theoretical gap with 
the transport gap measured in the intermediate coupling 
material YbBi2 yields good agreement, providing sup- 
port to the scenario presented for the field-induced gap 
closure. 
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